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Abstract 

An important ingredient of parton or string cascade models for ultrarela- 
tivistic heavy-ion reactions is a parton description of the baryon. Whereas 
previous models needed the concept of a diquark in an essential way, we have 
developed a new model using Dirac's approach of Poincaré-covariant many- 
body dynamics with constraints. 

In our model, the baryon is described as a dynamical set of three valence 
quarks and a fourth 'particle', the "junction", which carries the momentum 
fraction of the sea quarks as weh as aü of the glue. 

The model's parameters are the quark current masses, and one interaction 
strength, determined by the proton radius. Thus the model has no adjustable 
free parameters. Nevertheless, we obtain a remarkably good fit to the valence 
quark structure functions of the baryon. 
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Introduction. Among the methods employed in describing ultrarelativistic heavy ion 
reactions theoretically, various cascade models [0,11, (cf. also [§]) play an important róle. 
These models foUow the space-time development of the reaction dynamics in phase space; 
they are therefore essentially classical models. They simulate the nuclear reaction by a se- 
ries of elementary propagations, fusions and decays of partons or strings. When successful, 
cascade models not only reproduce the experimental data as asymptotic states of the re- 
action, but give some information on the dynamical, non-equihbrium development of the 
intermediate states, including a phase transition to a quark-gluon-plasma (QGP) — if such 
a phase transition should indeed occur. 

The basic objects in are strings; (confined) baryons or mesons are described by such 
strings sweeping (H-l)-dimensional periodic hypersurfaces in space-time ("yo-yo motions"). 
A baryon is thus modeled as a system of a quark and a diquark. Parton distributions are 
not included in this ansatz. 

In this letter we will present a new Poincaré-covariant dynamical model for a baryon, 
leading to yo-yo-hke motion of its constituents. Three valence quarks and a "junction" (rep- 
resenting sea quarks and all gluonic degrees of freedom), interact as classical point particles 
via a quasi-potential. This model, which is free of adjustable parameters, approximately 
reproduces the valence quark structure functions. The dependency of the structure 
functions will be discussed at the end of the letter. 

Constraint Dynamics. In constructing a Poincaré-covariant system of interacting 
point-particles, one has to face the consequences of the no-interaction theorem Q]. This 
theorem asserts that the only canonical Hamiltonian theory that is fully covariant under 
Poincaré transformations is that of a system of non-interacting particles. Dirac's method of 
introducing dynamical constraints [^] circumvents this problem. 

In this paper we use a set a (2A^) second-class constraints, written in the scheme of |^ 
(cf. also 0). The phase space is expanded to 8A^ dimensions : 



An evolution parameter s - without any direct physical meaning - is introduced for the 
parametrization of the phase space trajectories. A set of (2A^ — 1) Poincaré-covariant con- 
straints {$r \ r = 1 . . . {2N — 1)} and one connection between the evolution parameter 
s and the phase space variables are introduced which reduce the degrees of freedom of the 
system to a physical 6A^-dimensionaI hypersurface P' . We use the standard notation of 
writing A ^ B to denote that an equation holds only on F'; we also employ a summation 
convention over equal indices whereever applicable. 

The equations of motion are generated by a Poincaré-invariant 'Hamiltonian' H := X^^j 
via the Poisson brackets. 



The Lagrange parameters A-' are determined by the requirement that H preserve the con- 
straints: ^ ^ {$¿, $j} A-^ fti . Denoting the matrix of the Poisson brackets {$«, $j} 



l...N,fi = 0...3} 




(Ib) 



(la) 
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by Vij , the A-' are given by A-' ~ 
depends explicitly on s , we obtain 



A^ 



(V -^Y'^ ; since of all the constraints only 



-',2N d^2N 

ds 



(2) 



[It is straightforward to see that A^^ ~ ; thus H is independent of s] . Using Eqs. (1) 
and (H), one obtains the equations of motion for physically interpretable variables on the 
hypersurface T' . 

With constraint dynamics it is possible to construct a system consisting of two massive 
particles í7t.(i 2) , with equations of motion comparable to those of a string model It is im- 
portant that the period of this yo-yo motion is of the same order as that of the string model: 
(Tyoyo := 2M/KLund), bccause the period fixes the maximum time for one fragmentation [|T|. 
In order to achieve this, we use a quadratic quasi-potential (rather than the linear potential 
of the Lund model), which in turn increases the effective masses of the two particles: 



(m(i,2),g)' := (p(i,2))' = ("^(1,2))' - V 
= ("^(1,2)) - /í^ (g{i) - q{2) 
The dynamical system is defined by the constraints 

$1,2 := (p(i,2 

*3 := (g(i) - g(2)) • (P(l) +P(2) 



"^(1,2) j + /í (g{i) - g{2)^ 



3, ._ P(l)+P(2) ^ 

<I>4 .— ■ g(i) - s 

P{1) + P(2) 



(3a) 
(3b) 

(3c) 



The physical meaning of these constraints is the foUowing: $1^2 guarantee that the 
particles are on mass shell wherever the quasi-potential vanishes. Their effective masses 
increase with increasing distance . $3 ensures that in the center-of-momentum system 
(CMS) the particle distance is space-like , and $4 connects the evolution parameter s with 
the particle times in the CMS . 

The solutions of the equations of motion for this system are trigonometric functions (Fig. 
|l]). In the CMS and with vanishing particle masses, the period is T := ttM/{2k) , where 

:= (p(i) +P(2)) • Setting k = vr/4/íLund , one obtains the same period as in the Lund 
model. 

This fuUy Poincaré-covariant model can thus be thought of as simulating a meson of 
total mass M, composed of two quarks with current masses w^{i,2)- As in the Lund model, the 
quarks in the meson are confined, and this increases their effective masses. In comparison, 
the yo-yo string solution is only boost-invariant in one direction. Current quark masses are 
not included. 

Model of the Baryon. In the spirit of the above, a naive model for a baryon would be a 
system of three valence quarks, each carrying 1/3 of the baryon momentum. This, however, 
is too gross a simplification, since from deep inelastic scattering we know that about one half 
of the momentum is carried by sea quarks and gluonic degrees of freedom. In the model to 
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be described now, we will introduce a fourth (classical) particle, the 'junction' [1,0], which 
binds the valence quarks, and thus models these degrees of freedom (Cf. Fig. |). 

In keeping with the discussion in the previous section, we will again use a quadratic 
quasi-potential depending on the distances between the junction and the valence quarks 
[ again, constraints wiU ensure that these distances are space-hke in the CMS ] . This quasi- 

potential should affect the effective mass {p{,j)j of the junction, because pushing away a 
valence quark will increase the number of sea quarks and gluons, which are modeled by the 
junction. For simphcity, we will assume the valence quarks to always remain on mass shell. 
The above physical requirements are embodied in the following constraints: 

^»1,2,3 := (P{1,2,3)) - ("^(1,2,3)) (4a) 

i=i 

•^"5,6,7 := (^(1,2,3) - %J)) ■ P (4c) 
P 

*8 ■=^^-<l{J)-s , (4d) 

where P := + is the total 4-momentum of the system. 

The equations of motion for this system have no simple solution in closed form; but they 
are, of course, easily solved numerically. For the simulation of a proton, we use Mp = \P\ = 
0.938 GeV (proton mass), k = 0.5GeV/fm (leading to a proton radius of ~ 0.8 fm) and 



estimated values for the current masses of the quarks ||TT1: m(i) = m(2) = m(„) = 5MeV, 



m(3) = m((i) = lOMeV. 

Q"^ -independent Structure Functions. For fixed values of the variable x := Q"^ /2Mv, 
the substructure of nucleons is rather independent of the 4-momentum transfer (Q^). (M 
denotes the nucleon mass and v the transferred energy.) In the parton model this scaling 
behaviour is explained in terms of the presence of point-like charged constituents, called 
'partons'. Using the parton hypothesis, x can be interpreted as the fractional momentum 
of a parton in the infinite momentum frame. The distribution functions of this variable 
can be obtained from the data of deep inelastic e — p collision experiments. For comparison 
with our results, we shall use the parametrization of the valence quark distribution functions 
u.,{x),d^{x) given in [0] . 

To obtain valence quark structure functions in our model described by Eqs. (^), we first 
need the variable x in the context of our model. We define the longitudinal momentum 

fraction of particle i in terms of the ratio of the light-cone variables: X(j) := := pl_^p} 
. Note that this definition is fiavor-dependent because we are using fiavor-dependent quark 
masses (cf. end of the previous section). 

We have then sampled distributions of these momentum fractions by randomly choosing 
different initial conditions for the dynamical system of Eqs. (^; during the evolution of the 
system with one set of initial conditions, values of X(j) are sampled and averaged over. The 
results of these calculations are given in Fig. Rl and compared with the parametrization of 



m 



The rather good overall agreement of our result with the data seems rather remarkable, 
given the fact that ours is a purely classical particle model, and no free parameters were 
fitted in obtaining these results. 
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The tails of these distributions are definitely overestimated by our model. However, one 
should note that for x ~ 0.8 the number of events in the sample is vastly less than for x ~ 0.1 
(by about a factor of 10~^), so any error in the tail of the distribution hardly contributes to 
the momentum distribution of the proton. 

A fragmentation mechanism using the baryon model described here has been presented 
elsewhere |jl3|. 

Q'^-dependent Structure Functions. Experiments at HERA with high 4-momentum 
transfer show that the total momentum carried by the valence quarks decreases with in- 
creasing Q^, and that the number of gluons and sea quarks increases . In terms of 

our model, this would mean that the effective mass [P{j)] of the junction should increase 
with Q^, as it is the junction which models these degrees of freedom. However, since by 
construction the momentum transfer Q^ does not enter the dynamics of our model baryon 
at all, the only way to include this effect is by introducing an additional phenomenological 
mass term in the equations of motion for the junction. We thus replace the constraint <l>j 
in Eqs. (|) by 

*J := (P{J))'- ("^{J)(Q'))' + «:'E(?W ' (5) 

i=l 

with a phenomenological Q^-dependent term m(^j^{Q'^). Instead of simply fitting mi^j^{Q'^) 
to the data, we proceed as follows: The measured parton distribution functions for valence 
as well as for sea quarks and gluons, /fc(a;,Q^), are parametrized in From these we 

calculate integrated momentum fractions Xk{Q'^) ■= Jq xfk{x,Q'^)dx . Xgiue and Xsea are a 
measure of the amount of the corresponding 'mean field' carried by the J particle in our 
model and which we want to represent by m^j). However, since a certain part of this 'mean 
field' is generated dynamically by the system in the form of the effective mass squared of 

the junction, {p{j)j , we have to subtract this part, to avoid double counting. We are thus 



\ 2 

led to the ansatz 

m{j) 

where the proton mass Mp is needed for dimensional reasons. The factor a, finally, is 
determined by the requirement that rnf^j^^Q"^) vanish at the point where we connect with 
the Q^-iiidependent calculation: = 2 GeV^. Thus we take 'm(^j){Q'^) to be given by 



■^glue ~l~ -^sea ~^ ~ {•^dv + -^u^ , 



Mp 



2\ 



m{J)iQ ) := [Xg\ne{Q ) + a^sea {Q 

-^{¡^<iÁQ') + XuAQ'))]-Mp . (6) 

With this Q^-dependent term added, we again solve the equations of motion of the 
system and sample the distributions of the valence quark structure functions in the manner 
described above. The results are presented in Figs. || and |^, in comparison with [|14| . 



In conclusion, we have presented a truly Poincaré-covariant model for a baryon in terms 
of a completely classical particle description which, in spite of its simplicity, seems to de- 
scribe some of the features of the baryon surprisingly well. We believe this model (possibly 
with further refinements) to be useful in constructing models for ultrarelativistic heavy-ion 
reactions of the parton cascade type. We are currently in the process of developing such a 
code. 
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FIG. 1. Solution of the equations of motion in the CMS of the constraint dynamics model [Eqs. 
], compared to the Lund solution, in coordinate and in momentum space. í7i(i,2) = 0. 
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FIG. 2. Model for the baryon: a junction (J) is binding the valence-quarks (1,2,3) by a con- 
finement potential. The junction carries the average momentum of sea quarks and gluonic degrees 
of freedom. 



■o 

X 



3 

X 




FIG. 3. Sampled distributions of the valence quark longitudinal momentum fraction xuy{x) 
(up quark) and xd^^x) (down quark). The sohd and dashed hnes are the results obtained with the 
parametrization Íl! 
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FIG. 4. Sampled distributions of the valence quark longitudinal momentum fraction xu^{x) 
(up quark) and xdy{x) (down quark), at = (200 GeV)^. The soHd and dashed hnes are the 
results obtained with the parametrization 
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FIG. 5. Sampled integrated momentum fractions Xu,x¿ of the valence up quarks (□) and down 
quarks {o), as a function of Q^. The soüd and dashed Unes are the results obtained with the 



parametrization |14]. 
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